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ABSTRACT 

The failure of the XMM-Newton and Chandra X-ray telescopes to detect cooling gas in elliptical 
galaxies and clusters of galaxies has led many to adopt the position that the gas is not cooling at all and 
that heating by an active nucleus in the central E or cD galaxy is sufBcient to offset radiative cooling. In 
this paper we explore an idealized limiting example of this point of view in which hot, buoyant bubbles 
formed near the center return the inflowing, radiatively cooling gas to distant regions in the flow. We 
show that idealized steady state, centrally heated non-cooling flows can indeed be constructed. In 
addition, the emission-weighted temperature profiles in these circulating flows resemble those of normal 
cooling flows. However, these solutions are valid only (I) for a range of bubble parameters for which 
there is no independent justification, (2) for a limited spatial region in the cooling flow and (3) for a 
limited period of time after which cooling seems inevitable. Our exploration of non-cooling flows is set 
in the context of galaxy/group flows. 

Subject headings: galaxies: elliptical and lenticular, CD - galaxies: active - cooling flows - X-rays: 
galaxies - galaxies: clusters: general - X-rays: galaxies: clusters 



1. INTRODUCTION 

X-ray spectra of virialized hot gas associated with early 
type galaxies, galaxy groups and clusters indicate that the 
gas cools at less than a small fraction of the rate implied 
by their X-ray luminosities, (e.g. Peterson et al. 2001; 
Bohringer et al. 2002; Xu et al. 2002). There is no spec- 
tral evidence for gas at temperatures less than ~ T^riai/'^- 
As a result of these observations, it is often claimed that 
there is no cooling at all. Many authors have speculated 
that the gas is heated by an active nucleus (AGN) in the 
cluster-centered E or cD galaxy and that this heating can 
explain the absence of gas at T < T„iriai/3. 

The notion of AGN heating is also supported by the cen- 
trally disturbed gas visible in Chandra images of galaxy 
groups (e.g. Buote et al. 2003). Further evidence for 
hot bubbles on galaxy/group scales is provided by two (or 
limited multi) temperature spectral flts over an extended 
central region of the X-ray emitting gas in galaxy groups 
(Buote 2002; 2003; Tamura et al. 2003). The hot phase 
is typically 2 or 3 times hotter than the cooler phase Tc, 
but Tc(r) resembles the best- fitting temperature profiles 
in single phase spectral models. On the larger scales of 
galaxy clusters that contain powerful radio sources, huge 
bubbles of hot gas are apparently formed by strong cluster- 
centered radio galaxies and then buoyantly migrate out 
from the center of the fiow. (Fabian et al 2000; McNamara 
et al. 2000). Although the absence of apparent cooling and 
the evidence for X-ray surface brightness fiuctuations are 
shared by flows of all scales - galaxies or clusters - we 
are particularly interested here with galaxy/group flows 
because they may be less likely to have been disturbed 
by recent mergers and because the dominant stellar mass 
within the effective radius provides an independent mea- 



sure of the gravitational potential. 

The salient question that we address is: Can central 
heating greatly reduce the gas cooling rate M and also 
be consistent with the observed mean temperature of the 
hot gas? To explore possible answers to this question, we 
consider models in which bubbles of hot gas are heated to 
supervirial temperatures near the central AGN and then 
rise buoyantly to large radii. The bubbles move upstream 
in a normal cooling inflow in which gas slowly fiows toward 
the central potential minimum, losing energy by radiative 
losses. To avoid cooling to intermediate and low tempera- 
tures, the heated bubbles must carry outward all the mass 
arriving at the center of the cooling flow. The characteris- 
tic decline in the observed (single-phase) gas temperature 
toward the center of these flows, dT/dr > 0, is rather clear 
evidence that normal radiative cooling is driving an inflow 
of gas toward the center. In the circulation fiows we con- 
sider here, this cooling is assumed to be arrested by bubble 
formation associated with AGN heating. 

Flows in which cooling does not occur are expected to 
be in a steady state, at least on average. However, stellar 
mass loss and cooling or accretion inflow of hot gas from 
very large radii will cause the mass of gas in the circulation 
flow to increase slowly with time until some type of sud- 
den cooling occurs. However, it is possible that idealized, 
quasi-steady circulation flows may persist for sustained pe- 
riods without cooling, punctuated by intermittent cooling 
events. Remarkably, there is little or no observational ev- 
idence of ongoing or recent global cooling events in well- 
observed fiows. 

Can heated gas moving upstream in the cooling flow 
carry away all the mass arriving at the center and still be 
consistent with the low central temperatures and positive 
temperature gradients determined by (single-phase) X-ray 
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temperature observations? This is a question of some in- 
terest to us since our studies of AGN-hcatcd cooling flows 
using numerical gasdynamical simulations have been un- 
successful in agreeing with the central temperature profiles 
observed (Brighcnti & Mathews 2002; 2003). When AGN 
heating is sufficient to arrest central cooling, the gas cools 
at larger radii in the flow and the overall temperature gra- 
dient becomes negative throughout a large central region, 
unlike the observations. However, more analytic and ide- 
alized solutions without numerical diffusion, such as we 
consider here, may capture details of the flow that can- 
not be easily duplicated in gasdynamical simulations. To 
keep the apparent temperature from rising sharply toward 
the central AGN, it is necessary that the X-ray emission 
is not dominated by hot bubbles along the line of sight 
and that the inflowing gas is not heated too much by the 
counter-streaming bubbles. We show here that idealized 
steady state circulation flows can be constructed that arc 
consistent with observed temperature proflles, but only for 
bubbles in a restricted mass range and only for a limited 
spatial region and time. 

The discovery of X-ray holes or bubbles in rich clus- 
ters has resulted in a large number of numerical simula- 
tions. Most of these simulations focus on the evolutionary 
interaction of powerful radio jets and lobes from cluster- 
centered radio galaxies with the adjacent hot gas and do 
not quantitatively address the essential question whether 
or not strong radio jets can keep the gas from cooling to 
low temperatures and still preserve the density and tem- 
perature profiles observed. In the following section we 
briefiy review some of these numerical studies and con- 
trast them with the the heated flows we consider here. 

2. RECENT STUDIES OF AGN HEATING AND X-RAY HOLES 

The interaction of energetic supersonic jets with the hot 
gas has been simulated by many authors (e.g. Basson 
& Alexander 2003; Reynolds, Heinz & Bcgclman 2002). 
These numerical studies may accurately portray the initial 
active phases in which a high velocity jet from a central 
AGN flrst enters the hot gas. But supersonic jets create 
strong bow shocks that heat the ambient gas within a large 
cocoon. In most X-ray observations there is no c;onvinc;- 
ing evidence for extended shock-heated gas, although such 
phases may be relatively short-lived. 

In response to the apparent absence of gas at supcr- 
virial temperatures, many jet-lobe simulations have been 
designed that do not produce large masses of strongly 
shocked gas. These simulations typically involve an ad 
hoc slow release of energy (and mass) in fixed regions at 
some distance from the central AGN (e.g. Briiggen 2003a; 
Briiggen & Kaiser 2001, 2002; Briiggen et al. 2002; Chu- 
razov et al. 2001; Quihs, Bower & Balogh 2001). Briiggen 
and his collaborators have studied the consequences of in- 
jecting ultrahot gas with zero net velocity into fixed spher- 
ical regions of size ^ 1 kpc located ^ 10 kpc from the cen- 
tral AGN in galaxy clusters. Evidently, this type of energy 
release is an approximation to a stationary working sur- 
face of a jet. The hot gas atmosphere surrounding the 
AGN is initially in hydrostatic equilibrium with outward 
increasing entropy. The luminosities of the energy injec- 
tion, L ~ 10**^ — 10**'^ erg s~^, are usually chosen for each 
simulation to be near the maximum possible L that can 



still produce buoyant plumes and bubbles; larger L would 
evacuate gas throughout the entire cluster core in times 
shorter than the buoyancy rise time (e.g. Briiggen 2003a). 
A common feature of these solutions is that cooler, low 
entropy ambient gas follows upward in a transient "en- 
trained" flow behind the rising hot gas. This transient, 
approximately axisymmetric convective circulation occurs 
as the rising gas is replaced by downflowing ambient gas 
surrounding the jet axis that turns and flows upward be- 
hind the heated gas. After ~ 100 Myrs, when the energy 
injection is turned off, some of this "uplifted" low entropy 
gas typically begins to fall back toward the center, perhaps 
re-energizing the central AGN at a later time with a feed- 
back mechanism. While these papers are useful contribu- 
tions to our understanding of the interactions of powerful 
radio jets with the hot gas, the authors often note that the 
cooling flow M may be reduced if energy is released in this 
manner with duty cycle ^ 0.5. However, the mass cooling 
rate M{t) near the flow center is not explicitly calculated 
so it is difficult to link these simulations to the problem of 
the low upper limits on M set by XMM observations. 

By contrast, our numerical simulations study feedback 
heating and M reduction in genuine cooling flows over 
times comparable to the Hubble time (Brighenti & Math- 
ews 2002, 2003). The original gas is heated by the central 
AGN and no new mass is introduced. M{t) at the flow cen- 
ter and elsewhere is monitored throughout the entire cal- 
culation. In our experience, when these cooling flows are 
sufficiently heated to reduce the central cooling rate M, as 
XMM observations require, the azimuthally-averaged tem- 
perature and density profiles typically deviate from those 
observed. When cooling is suppressed by heating near the 
center, we find that it occurs at the same rate but at larger 
radii in the flow. In flows with off-center jet-like heating, 
the gas at large angles relative to the heating axis easily 
flows toward the center and eventually cools to low tem- 
peratures. 

The flow heating problem on galaxy/group scales that 
we consider here is very different from the hciating of hot 
gas in richer galaxy clusters by powerful, linearly oriented 
radio jets and lobes. Many flows with undetected M do not 
contain powerful radio jets. For example, when observed 
with XMM- Newton, the galaxy/group NGC 5044 shows no 
RGS or EPIC spectral evidence of gas cooling to temper- 
atures below about 0.7 keV (Tamura et al. 2003; Buote et 
al. 2003). The (single phase) gas temperature profile drops 
from about 1.2 keV at 30 kpc down to 0.7 keV at the origin, 
similar to many other E galaxies and qualitatively similar 
to T(r) in rich clusters. Also, like most elliptical galax- 
ies, NGC 5044 has an unexceptional and unresolved weak 
central mJy GHz point source (NRAO/VLA Sky Survey; 
Condon et al. 1998). However, the Chandra image of 
NGC 5044 (Fig. 1 of Buote et al. 2003) shows that the 
hot gas within at least ~ 14 kpc is highly disturbed in 
an azimuthally random and chaotic manner. In addition, 
temperature fluctuations in the range 0.7 ^ kT < 1.4 keV 
are observed out to ~ 30 kpc (Buote et al. 2003). The 
most important morphological attribute of these thermal 
and density irregularities is that they are symmetric about 
the center of NGC 5044, with no evidence of jet-like ori- 
entation. If AGN heating is the agency that keeps the gas 
from cooling, the heating in NGC 5044 and other simi- 
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lar E galaxies must be more radially symmetric than the 
jet-heating considered in the radio galaxy simulations dis- 
cussed above. Alternatively, because of the short dynami- 
cal times in galaxy scale flows, the jet axis in would have to 
be randomly reoriented every ^ 10^ years to maintain the 
appearance of symmetric heating. Finally, the outward 
mass flux M is likely to be maximized in spherically sym- 
metric bubble outflows that move away from the central 
AGN in every direction. 

Ruszkowski & Begelman (2002) and Briiggen (2003b) 
have studied non-cooling, static models of cluster gas at- 
mospheres. The gas is heated at large radii by thermal 
conduction and at small radii by a cloud of tiny expand- 
ing bubbles, a heating mechanism proposed by Begelman 
(2001). The heating is spherically symmetric, similar to 
the flows we discuss here. These cluster-scale models are 
almost uniquely successful in preserving the positive tem- 
perature gradient observed. However, this combination of 
central heating and conduction does not work for flows 
on galaxy/group scales; even with the maximum (Spitzer) 
conductivity it is impossible to get both small M and 
acceptable temperature gradients (Brighenti & Mathews 
2003). In the heating scheme proposed by Ruszkowski 
& Begelman a large number of very small "effervescent" 
bubbles are produced near the central AGN and buoyantly 
flow into the hot gas throughout a large volume, acquiring 
a steady state radial distribution. The bubbles do not mix 
thermally with the ambient gas but expand adiabatically, 
heating the ambient flow by PdV work as they expand. 
While this bubble model has some features in common 
with that considered here, we find that very small bubbles 
are inconsistent with steady state circulating flows. 

In the following discussion we investigate the possibil- 
ity that hot gas in galaxy/groups similar to NGC 4472 or 
NGC 5044 is heated near the center of the flow and con- 
verted into hot bubbles. The rising bubbles are required 
to carry mass outward at the same rate that mass arrives 
in the core by normal radiative losses in the interbubble 
flow. The interbubble flow receives both momentum and 
energy from the rising bubbles. To optimize the outward 
rate of mass transfer, the heating is assumed to be sym- 
metric about the central AGN, not conflned to a single 
jet. 

3. CIRCULATION FLOWS 

We explore here the conditions required to estab- 
lish oppositely moving, two phase, steady state flows in 
galaxy/groups in which the gas does not cool to low tem- 
peratures. A cool phase near the virial temperature flows 
inward as in a normal cooling flow. At some small radius 
rh we assume this flow is heated and transformed into hot 
bubbles that buoyantly float upstream in the cooling flow, 
returning the mass flux received in the incoming flow to 
some large circulation radius Vc where the bubble entropy 
matches that of the local cooling flow gas. Gas circulates in 
the gravitational fleld like a thermodynamic engine. The 
mechanism that produces the bubbles is not specified and 
for simplicity the physical condition of the gas within or 
beyond Vc are not considered. Clearly, an idealized steady 
flow model like this cannot be strictly correct since the 
mass of circulating gas is expected to gradually increase 
due to stellar mass loss and inflow of gas from r > Tc- As 



the mass of circulating gas slowly increases, the increas- 
ing gas density should eventually lead to extensive cooling. 
Nevertheless, a quasi-steady circulation model could be a 
reasonable representation for a limited time between cool- 
ing episodes. Unlike the numerical simulations described 
earlier, no new mass is injected into the flow. 

We also assume that the hot bubbles retain their co- 
herency as they move through ~ 10 — 15 bubble diameters. 
This degree of coherence is in fact observed in the Perseus 
Cluster (Fabian et al. 2000), assuming that the distant 
"ghost cavities" were formed closer to the central AGN. 
The hot bubbles simulated by Churazov et al. (2001) float 
upward in the cooling flow atmosphere and come to rest 
at some large radius where the ambient entropy is equal to 
that within the bubble. However, we have found bubble 
evolution to be a diflicult problem for numerical gasdy- 
namics and we remain sceptical of our assumption of bub- 
ble coherency. The bubbles in Perseus may be stabilized 
by internal magnetic fields. It would be relatively easy to 
modify our circulation model to include bubble fragmenta- 
tion according to some ad hoc prescription, but the effect 
of bubble fragmentation can easily be imagined by com- 
paring similar circulation flows having different constant 
bubble masses. 

In our circulating fiows the hot bubbles do not thermally 
mix (by thermal conduction for example) with the incom- 
ing flow before reaching the circulation radius Tc where 
the bubble and inflow entropies are approximately equal. 
Thermal mixing is another effect that could be included 
in more complicated circulation flows, but our assumption 
of no thermal mixing is a limiting case of interest. If pre- 
mature bubble-flow thermalization occurs, the inflowing 
gas will be heated and the apparent temperature profile is 
likely to become negative, unlike the (single phase) tem- 
perature profiles observed. Unintended thermal mixing 
may have occurred in our numerical simulations of heated 
cooling flows because of the limited numerical resolution 
(Brighenti & Mathews 2002; 2003). When the bubbles and 
interbubble flow are viewed along the same line of sight, 
an apparent positive dT/ dr can be maintained if the filling 
factor /h of the rapidly moving bubbles is relatively small 
and if the bubbles do not significantly heat the ambient 
gas. We show that the constraint /b < 1 also limits the 
heating that bubbles can receive at the heating radius rh 
and therefore also the radial extent Tc of the circulation. 

Finally, to keep our circulation models as simple as pos- 
sible, we have resisted a temptation to include an addi- 
tional pressure component within the bubbles due to (rel- 
ativistic) cosmic rays or magnetic fields. An additional 
nonthermal pressure could help increase the outer radius 

of the circulation, but may transport k;ss mass outward. 
Nevertheless, many of the same constraints we derive here 
for purely thermal circulations also apply to bubbles with 
nonthermal buoyancy. 

Our approximate, idealized circulation fiows are not in- 
tended to reproduce the observed properties of any partic- 
ular group or cluster although we use parameters for the 
massive elliptical galaxy NGC 4472 whenever appropriate. 
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4. INFLOW EQUATIONS 

The steady state equation of continuity for the inflowing 



gas 



M = Aur fpu 



(1) 



contains the fraction of the volume /(r) fiUed by the in- 
flowing gas. The bubble filling factor is /b = 1 — /. Both 

M and ti are negative. The equation of motion and the 
conservation of thermal energy in the cooling inflow are 



du 

pu— = — 
ar 



d{p0) 
dr 
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and 



3^6* ^ dp p^k 

pu-— = OU- TT- 



2dr 



dr 
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Here we define 6 = kT/prUp where p is the molecular 
weight so the gas pressure is P = p9. The radiative cool- 
ing coefficient A(T, z) erg cm"^ s~^ is based on the results 
of Sutherland & Dopita (1993). The effective gravity ge 
includes the momentum imparted to the inflowing gas by 
bubbles moving in the opposite direction (see §5 below). 

These equations can be written as two coupled differen- 
tial equations, 



du _ u [5a9 - '3ge + {2rpK/mpU)] 
dr ~ r (3u2 - 561) 



and 



d9 _ 2 [ue{ger - u^a) + {rp{A/ml){e - u^)] 



dr 



(4) 



(5) 



r (3u2 - 50) 

and p{r) is found from Equation (1). In deriving these 
equations, dp/dr has been eliminated by differentiating 
Equation (1) and the variation of the filling factor is in- 
cluded in the variable 

o , dlogf 



dlogr 



(6) 



The gravitational acceleration g is based on the massive 
E galaxy NGC 4472 which is embedded in a galaxy group 
having an NEW dark matter halo. The stellar mass has 
a de Vaucoulcurs profile (total mass: M^t = 7.26 x 10^^ 
Mq; effective radius: Re = 1.733' = 8.57 kpc) with a core 
P*,core(?") = P*,deV {rb){r /rb)~°-^° within the break radius 
Vh = 2.41" = 200 pc (Gebhardt 1996; Faber et al. 1997). 
The NEW halo has total mass M/^ = 4 x lO^-"' Mq. For an 
excellent fit to the gravitational force in this stellar-dark 
matter potential we use 
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cm s 



(7) 



with a = 1.0233 x 10^^^ ^ ^ q.I, b = 4.0378 x 10"^, 
q = -0.8354 and s = 3. 

The hot gas density observed in NGC 4472 is found 
by combining data from Einstein (Trinchieri, Fabbiano, & 
Canizarcs 1986) and ROSAT (Irwin & Sarazin 1996). Eor 
the inner region we have Abel-inverted Chandra surface 
brightness data from Loewenstein et al. (2001) and nor- 
malized them to previous observations. These data can be 
accurately fit with 
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n{i) 



+ [r/i?er(i)]pW} 



(8) 



(9) 



where 

n{i) = (0.117,0.00735,-0.000493,0.142) 
r{i)/Re = (0.107, 0.95, 10., 0.04) 
p{i) = (2.0,1.14,1.19,3.0). 

The gas temperature in NGC 4472 

Tir) = 1.55 x 10^ IM^±1:L 
rkpc + 5.98 

is a fit to the combined ROSAT HRI and PSPC observa- 
tions of Irwin & Sarazin (1996) and ROSAT PSPC obser- 
vations of Buote (2000). 

4.1. Flow Without Bubbles 

The cooling flow equations are solved indirectly by re- 
quiring that the gas density match the observed density in 
NGC 4472 at the end points of the circulation, r/j and rc- 

To illustrate a typical solution of this two-point boundary 
value problem, we consider a circulation flow between radii 
rh = 1 and Tc = 10 kpc. The instantaneous rate that gas 
cools in this annular region can be determined from 



M 



p'^A 2/xmj 



Anr'^dr w 0.7 Mq yr'^ (10) 



using the observed T(r) and p{r) proflles and evaluating 
A with an assumed abundance 0.5 solar. 

Figure 1 shows a typical solution of Equations (1), (4) 
and (5) for a cooling flow without bubbles (/ = 1, a = 2, 
and ge = g) between ri and r2 assuming M = 0.7 Mq 
yr^^. When the flow equations are solved iteratively in 
the manner described, both the temperature proflle and 
its end-point values are completely determined by p{rh) 
and p{rc)- The modest disagreements apparent in Fig- 
ure 1 between the computed temperature and density and 
the proflles observed in NGC 4472 are expected because 
of the approximations we have made: steady state flow, 
suppression of mass loss from stars, no gas extending be- 
yond rc, etc. When the flow is subsonic [u^ <C 6) the 
temperature gradient dO/dr (Equation 5) has the same 
sign as uge + pK/rn^. Since w < 0, it is clear that the 
radiative term is responsible for the positive temperature 
gradients observed. Although the flow temperature gradi- 
ent in Figure 1 is flat, dT/dr « 0, this solution is entirely 
adequate for our purposes here: to determine the influ- 
ence of counter-moving bubbles on the observed thermal 
properties of the flow. We are primarily interested in dif- 
ferential changes in the temperature profile when bubbles 
are introduced. 

5. BUBBLE DYNAMICS 

A circulation flow without cooling below T{rh) can be 
maintained if the inflowing gas is heated at radius rn, 
forming bubbles of hot gas that rise buoyantly upstream 
in the cooling flow. At Vh the bubbles have temperature 
Tb{rh) = hT{rh) where T{rh) is the local cooling flow tem- 
perature and /i > 1 is a dimensionless heating factor. To 
avoid needless complications, all bubbles are assumed to 
have the same mass ruh = p6(4/3)7rr^ for each solution; the 
consequences of bubble breakup can be inferred from so- 
lutions with different mt,. As they move outwards, the hot 
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bubbles evolve adiabatically (7 = 5/3) in pressure equilib- 
rium with the ambient gas. The density inside the bubbles 
is 

h \e~h) 



Pb 



(11) 



where ph = piru) and Oh = 0{rh). 

The equation of motion of the bubbles is dominated by 
the drag force so the terminal velocity is reached almost 

immediately after they are formed. When moving at the 
terminal velocity, the forces on the bubble balance: 

dP Aru 

mtg + + 5pf{u, - ufnrl = (12) 

where Uf, is the bubble velocity. The external pressure 
force on the bubble is the pressure gradient across the bub- 
ble dP/dr times the mean chord through a sphere 4rb/3 
and the bubble area Trr^. The final term in Equation (12) 
is the approximate drag force on the bubble where p is the 
density in the ambient cooling flow and (5 is a dimension- 
less coefficient of order unity that depends on the shape of 
the bubble and other imponderables; we assume 5 = 0.5. 
Since the cooling flow is highly subsonic, hydrostatic equi- 
librium is a good approximation, dP/dr « —pg, and the 
bubble velocity follows from Archimedes' principle, 



Ub = U + 



6(3 \pb p 



1/2 



where 



rub rub \4:npbJ 



2/3 



(13) 



(14) 



The collective drag interactions of a multitude of rising 
bubbles exerts an outward force on the cooling flow so the 
effective gravity becomes 

ge = g- Siub - uf{l - f)pbp. (15) 

Typically, ge differs from 5 by < 5 percent. 

The energy transmitted from the bubbles to the sur- 
rounding gas is of particular interest for the cooling flow 
problem since the observed temperatures are always lowest 
at the center (as in Equation 9). This observation indicates 
that most of the X-ray emitting gas near the flow center 
is not being strongly heated even when heating is required 
to account for the absence of further cooling. 

The nature of the energy exchange between the bub- 
bles and the flow is not well understood. Several authors 
(Briiggcn & Kaiser 2002; Churazov et al. 2002 Bcgelman 
2001; Briiggen 2003b) have suggested that cooling flows 
are heated by the PdV work expended as bubbles expand 
during their ascent in the flow. While some PdV heat- 
ing by bubbles may occur during the transient adjustment 
when bubbles first move into a region of the flow devoid 
of bubbles, it is less obvious that strong PdV heating is 
expected in steady state flows where the filling factors of 
the bubbles and flow remain constant at every radius, i.e. 
dV = 0. In steady flow as expanding bubbles rise subson- 
ically upward into a spherical shell of interbubble gas, the 
volume of interbubble gas in that shell has already been re- 
duced to accommodate the expanded bubbles; this volume 
adjustment is contained in the filling factor /(r) for that 
shell. Since the bubbles move subsonically, the flow in the 
interbubble gas as it is rearranged at constant volume in 



the shell is also subsonic and therefore largely incompress- 
ible, i.e., dV = 0. PdV heating is likely to depend on the 
filling factor in some complicated manner that we do not 
investigate here. In spite of these reservations, we consider 
the possibility that PdV heating may be important. 

In steady state a single rising bubble of increasing vol- 
ume Vb = rub/pb delivers PdV power at a rate 

Epdv = Pb-^ = Pb{Ub - ")-^ 

TUb r dpb , -.rribpg _i 

= -Pbiub-u) — = iub-u) ergs 

PbT Pb dr ^pb 

where we assume P = Pb cc pj and dP/dr = —pg. Multi- 
plying by the space density of bubbles, Nb = — f)pbl'>Tib) 
wc find the power delivered to the inflowing gas by all bub- 
bles. 



^^{'^-f){ub-u)pg 
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where epdv is the uncertain efficiency that the expanding 
bubbles heat the gas locally. 

Regardless of the uncertain relevance of PdV heating, 
the cooling inflow is certainly heated by the drag interac- 
tion with the counterstreaming bubbles. The drag-related 
power generated by a single bubble is the product of the 
drag force and the relative velocity Fd{ub — u) = 5p{ub — 
7i)^7rr^. The local density of bubbles is A^5 = pb{l — ,f)/mb 
so the collective power density produced by all bubbles is 
Fd{ub — u)Nb. The bubbles could heat the surrounding 
gas by producing (weak) shocks or turbulence that dissi- 
pates locally. Non-local heating could result if the energy 
received from bubbles was in the form of sound waves, 
gravity waves or long-lived turbulence. Since the fraction 
of this energy that goes into local heating is unknown, we 
include a heating efficiency factor Cd in the power den- 
sity. The heating efficiency of rising bubbles is discussed 
by Churazov et al. (2002). Bubble-flow drag interactions 
heat the gas at a rate 
3ed 

Sd = -. < 

in 

Bubble heating can be included in the flow equations by 
replacing the radiative cooling coefficient A in Equations 
(4) and (5) with A -|- Apdy + Ad where 



erg cm s "'^ 



-S{1 — f)p{ub — u)^ erg cm ^ s ^. 



Apdv 

and 



3 ^giub-u) 
o p 



erg cm^ s ^ 



Ad = --edS{l- f)—^- ergcm'^s^ 

4 p n 



(16) 
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5.1. Steady State Circulation Flows 

The essential and necessary condition for steady state 
circulation requires that the rate of cooling mass inflow 
M is balanced at every radius by a proportional upstream 
mass flow M carried by the bubbles, 

M = Airr^fpu = -AT^r'^{l - f )pbUb (18) 

so that 



/ = 



PbUb 



(19) 



PbUb - pu 

The luminosity expended by AGN heating at radius rh, 

Zpnip 



6 



M 



10^ Ky I Mq yr-i 



erg s 



can be produced by accretion onto the central black hole 

with mass flow Albh ~ Lagn/rjc^ where ry w 0.1 is the ef- 
ficiency that the accreted mass converts to energy. Com- 
bining these expressions, 



T 



M 



W KJ \ Mq yr-1 



Since Mbh "C M, only a negligible flow into the hole can 
provide the necessary heating at rh so the mass flux in 
both radial directions in the flow beyond Vh must exactly 
balance at every radius (Equation 18). 

Counterstreaming flow solutions are found by iteratively 
solving Equations (1), (4), (5), (11), (13) and (19). We be- 
gin with a bubble- free flow as in Figure 1. then calculate 
the dynamics of the bubbles required to carry the same 
M . The total volume of the bubbles at each radius pro- 
vides an estimate of the filling factor /(r), which can then 
used to update the (iterative) solution for the inflowing 
gas. The iteration of the two point boundary problem for 
the cooling inflow (Equations 4 and 5) is embedded within 
the global iterative scheme involving the bubble dynamics. 

For simplicity we fix the inner radius rh of the circu- 
lation for all calculations, but the outer radius rc is free 
to vary as required by the solution. With each bubble- 
flow iteration, the radius of the circulation flow Tc is de- 
termined as that radius where the entropy of the rising 
bubbles equals that in the cooling inflow; the cooling flow 
density is then rcnormalizcd to the local density of NGC 
4472 at Tc. At r — Tc the density and temperature inside 
the bubbles naturally become the same as those in the 
local flow, and the bubbles merge with the cooling flow. 
The outward velocity of the bubbles at Tc and the inward 
velocity of the cooling flow gas are both very small and 
highly subsonic at but never exactly zero; we assume 
that this small differential velocity is dissipated without 
appreciably altering the local entropy. As before, the tem- 
perature profile of the interbubble flow and its end-point 
values T{rh) and T{rc) are completely determined by the 
requirement that the gas density match that observed in 
NGC 4472 at Vh and r^. 

For a rigorously consistent solution the global mass flow 
rate M should be redetermined from Equation (10) (using 
the computed flow) for each value of rc as it changes dur- 
ing the iterative calculation. However, we have found that 
this requirement is not essential for overall convergence 
and for simplicity we have decided to relax this additional 
constraint, keeping both r/, = 1 kpc and M = 0.7 Mq yr~^ 
constant for all solutions. For given M and rh, the impor- 
tant parameters for each circulation flow are the heating 
factor h, the bubble mass mb and the heating efficiencies 
Cpdy and €4 of the bubbles. 

6. RESTRICTIONS ON FLOW PARAMETERS 

Not every set of flow parameters h, mb, e^, Cpdv, M 
and rh - results in a solution that converges. Successful 
solutions must satisfy several additional conditions: (1) 
the bubble radius must not exceed the local radius in the 
flow, (2) the bubbles must not fill the entire volume and 



choke off the inflowing gas, and (3) the mass of the bubbles 
must be sufficient to carry the outward mass flux M. 

The first obvious restriction is most acute near the heat- 
ing radius rn where rb/r is largest. The requirement that 
the size of the bubbles not exceed r/j is 



rh> rb = 0.44 



rub 



1/3 



kpc. 



^(106 Mq) {n,{rh)/0.1 cm-' 
The second condition is that the filling factor for the 
bubbles fb = ^ — f cannot exceed unity near rh where 
the bubble density and flow congestion is highest. If 
tbd = '^rb/[ub{rh) —u{rh)\ is the time for bubbles to move 
one bubble diameter at r = rh, tb = nib/M is the time 
interval between bubble formation, Vb = (4/3)7rr^ is the 
bubble volume, and Vh « STrr^rf, is the volume of an an- 
nular shell of thickness 2rb <C rh, then we require that 



(20) 



^ 17 (21) 



This condition can be written 

M h 

{Mqjv-^) MbhT','\l^^Mh/Q.lc 

where Mhh = [w6(?'h)— M(?'fe)]/(56'ft/3)^/^ is the Mach num- 
ber of the relative flow and Uh = {ph/mp){2 + fi)/{5iJ,). 
Condition (21) sets an upper bound on the amount of heat- 
ing h that can occur at rh- 

Finally, from Equation (18) the mass flux carried by the 
bubbles is M = 47rr^(l — f)pbUb- Near the heating radius, 
r ^ rh, where pb < p and Ub > \u\, the bubble velocity 
Ub ~ {g/Sppby/^. Therefore, 

1/3 , 5/6 . g ^ 1/2 



M«47rr^[l-/(r^) 



\Tr6 



ra 



1/6 



• ^^^^ 

where pb ~ ph/h. For a given cooling flow where M, rh, 
ph, f and g are all specified or constrained, the bubble 
mass jrib must be large enough to carry the mass flow M 
dictated by the cooling inflow. Although the dependence 
on mb is weak, mb and h are the only flexible parameters. 



As mb decreases, the central heating factor h cx 



1/5 



must 



also decrease for fixed M. If bubbles of very small mass 
rrib are formed at rh or are produced later on by the frag- 
mentation of larger bubbles, condition (22) can cause the 
maximum allowed heating h to fall below unity for which 
no circulation flows are possible. 

7. RESULTS 

7.1. Flows Without Bubble- Flow Energy Exchange 

Results for a sample of steady circulation flows are il- 
lustrated in Figure 2. The heating efficiencies and Cpdv 
that the flow is heated by the bubbles are assumed to be 
zero and the mass flow rate M = 0.7 M© yr~^ and heat- 
ing radius rh = 1 kpc are constant for all flows. Two 
circulation flows are shown in Figure 2 for each assumed 
bubble mass mb, one with a typical value h and one near 
the largest h = hmax for which solutions are physically 
possible or appropriate. Each solution in Figure 2 shows 
the cooling inflow variables n^ir), T{r) and u{r) with solid 
curves. In the density panels the solid curves begin and 
end on the observed density profile for NGC 4472, which 
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is shown with dotted hncs in the top row of panels. When 
this flow reaches rh — I kpc, the gas is heated by a fac- 
tor h and the density and temperature change abruptly 
to bubble values, pbirh) = p{rh)/h and Th{rh) = hT{rh) 
respectively. Then the bubbles rise in the flow along the 
dashed lines until they eventually intersect the cooling in- 
flow at some large radius rc- The entropy of the bubbles 
and flow are identical at Tc and the bubbles provide just 
enough gas to maintain the inflowing cooling flow from 
that radius. The third row of panels (c, g and k) shows 
the negative cooling inflow with solid lines and the much 
faster outward motion of the rising bubbles, «(, > 0, with 
dashed lines. The variation of the bubble radius r6(r) is 
shown in the bottom row of panels (d, h and t) with dash- 
dotted lines. The dashed lines in these bottom panels show 
the variation of the inflow volume fllling factor /(r). As h 
approaches the highest possible value for each bubble mass 
mf,, h — > hmax, the bubbles nearly fill the entire available 
volume at where / — > 0, as explained by Condition (21). 

For the most massive bubbles considered, mb = 10^ Mq , 
solutions for two heating parameters are shown, h = 3 
and 6 for which the circulation radii are rc = 10 and 
25 kpc respectively. The density, temperature and veloc- 
ity profiles for the cooling inflows for these two solutions 
are almost identical. We have not considered larger val- 
ues of h because the bubble velocity Ub for the h = 6 
solution approaches the sound speed in the cooling flow 
(561/3)1/2 ^ (r/io7 K) km s'^ Larger h would re- 
sult in supersonic bubble velocities that would drive shocks 
into the cooling gas, dramatically increasing its tempera- 
ture. Also the bubble size at r/, for the h = 6 solution, 
fb(jh) ^ 1 kpc is comparable to r^, again arguing against 
larger values of h. 

The central column of Figure 2 illustrates two circula- 
tion flows for bubbles of mass mft = 10^ Mq for /i = 3 and 
6.5. The circulation with h = 6.5 extends out to Tc = 18 
kpc, but the flow filling factor / is less than about 0.1 at 
rfi so there are no circulations with h > 6.5 ~ hmax (Con- 
dition 21). In order for this marginal flow to support mass 
flux M, the relative flow Ub—u must increase as r ^ r^i and 
the drag increases accordingly. The bubble velocity Ub is 
seen to decline toward rh- As a result of the large drag, the 
effective gravity ge (Equation 15) is directed outward near 
Th = I kpc where the flow would be Rayleigh- Taylor un- 
stable. This region of instability becomes larger for higher 
values of h. Because of the approximations we have made, 
the density in the flow region < r < Vc is somewhat 
higher than the density observed in NGC 4472 (dotted 
lines). However, if this flow were observed, the apparent 
density rig^app would be lowered since the flow only occu- 
pies a fraction / of the volume, i.e., Ue^appif) ~ f^^'^neir). 

The first column in Figure 2 shows the circulation pat- 
tern for bubbles of mass mb = 3 x 10* Mq for two flows 
with h — 3 and 3.6. The flow fllling factor / is drop- 
ping rapidly as r — + r?j for the h = 3.6 solution so this is 
close to the maximum possible heating (and rc) for these 
smaller bubbles. We were unable to achieve fully con- 
vergent flow solutions with even smaller bubbles of mass 
rUb ^ Mq. While our search for such solutions was 
not exhaustive, it is apparent from the trend in panels k 
— > 5 — > c in Figure 2 that the maximum values ofub, h and 
Tc all decrease with rub, as expected from condition (22). 



Bubbles of mass nib = 3 x 10"* M© are approaching the 
minimum mb necessary to carry the cooling mass flow M 
for any heating h > 1. Non-circulating flows with bubbles 
of lower mass must ultimately cool radiatively to tempera- 
tures much lower than T(r/j), in disagreement with XMM 
observations. 

One of the motivations for our exploration of non- 
cooling circulation flows is to estimate the influence of the 
bubbles on the apparent gas temperature which typically 
has a minimum near the origin. To simulate the apparent 
temperature that would be inferred when viewing the cool- 
ing inflow and the hot rising bubbles along the same line of 
sight, we calculated the local emission- weighted tempera- 
ture 

p^Tf + plTb{l-f) 



(T) 



P'f + plil-f) 



(23) 



which is plotted as dotted lines for each of the flows in pan- 
els 6, / and ) of Figure 2. For most circulation flows (T) is 
nearly indistinguishable from the temperature of the cool- 
ing flow component. Consequently, a single phase inter- 
pretation of the temperature profile would show a mildly 
positive gradient, dT/dr > 0, similar to those observed 
and consistent with traditional cooling flows. This is a de- 
sirable feature of the circulation model. Only for the least 
massive bubbles considered, Wb = 3 x 10* Mq, does (T) 
rise slightly as r ^ r^t, but the differential change in the 
temperature proflle caused by the hot bubbles is small. 
The apparent density papp = ((p^))^^^ and entropy pro- 
files are also similar to those of the cooling flow component 
alone. We stress that these results differ from normal con- 
vection where (1) there is a single temperature and density 
at every radius (apart from fluctuations), (2) the tempera- 
ture gradient is determined by the gravitational potential, 
dT/dr = —{2pmp/5k){d^/dr), and therefore has a sign 
opposite to that observed, and (3) the entropy is constant. 

7.2. Flows With Bubble-Flow Energy Exchange 
The ratio of PdV heating to drag heating is 
ipdv _ 4 Cpdv grb _ ^ ^ (rfc/lOOpc) 
id 5 Cd 5{ub-uY 



1.0- 



'^%c K/lOO km/sec)2 



-0.85 
kpc 



cm 



where Ub <C \u\ and we use g ~ 4.4:7 x 10 
s~^ which is appropriate for 1 <> rkpc ^ 100 in NGC 4472. 
The result ipdv/id ^ 1 follows immediately from Equation 
(13) provided Ub <C |m|. In view of the uncertain physical 
nature of these two comparable bubble-flow heating mech- 
anisms, we consider heating due to bubble expansion and 
drag separately. 

Figure 3 illustrates four representative cooling flow mod- 
els in which the work done by the drag force is assumed to 
heat the local cooling flow gas with various efficiencies 
(and Cpdv = 0). As before M = 0.7 M© yr~i and rn = I 
kpc for all flows. The four flows also have the same bub- 
ble mass mb = 10^ Mq and heating factor h = 6, but the 
bubble drag heating efficiency has values = 0, 0.1, 0.3 
and 0.7. Flows with progressively higher T, Tb and n near 
the heating radius r^ correspond to increasing efficiencies. 
By contrast, the circulation radii, velocities, filling factors 
and bubble radii are insensitive to the heating efficiency. 

However, the flow with the largest heating eSiciency, 
Bd = 0.7 is less satisfactory because flow temperature and 
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particularly the apparent temperature (T) develop appre- 
ciable negative radial gradients near rh- Because of this 
rise in the flow temperatures, bubbles heated with the 
same h rise to higher temperatures. In addition, when 
= 0.7 the density gradient develops an unstable posi- 
tive density gradient within about 2 kpc. From these and 
other circulation flows with nonzero we conclude that 
the apparent temperature has no undesirable radial gra- 
dients as long as ^ 0.5; this is a reasonable constraint 
in view of the uncertain nature of the bubble-flow heating 
process. 

Figure 4 shows three flows assumed to be heated only by 
bubble expansion with efficiencies Cpdv = 0.1 0.5 and 0.7 
(and 6(1 = 0). As before, the apparent emission- weighted 
temperature profiles (dotted lines in Fig. 4b) become 
strongly negative unless Cpdv ^0.5. 

8. FINAL REMARKS 

Our primary objective in constructing these models is to 
explore the possibility that cooling flow gas can be heated 
at some small radius in such a manner that (1) no gas cools 
to very low temperatures and (2) the apparent tempera- 
ture gradient of the cooling flow does not become negative 
in the central regions. We have demonstrated that flows 
satisfying these two criteria are indeed possible. With ap- 
propriately chosen parameters, buoyant bubbles can carry 
the cooling flow mass M away from the core with velocities 
sufficiently large that the apparent (emission-weighted) 
temperature gradient d(T)/dr is virtually imaffcctcd by 
the hot gas inside the bubbles. In addition, the local in- 
terbubble gas can be heated by the bubbles with up to 25 
- 50 percent efficiency without adversely disturbing the 
temperature gradient. This positive assessment of ideal- 
ized circulation flows contrasts with the less satisfactory 
results of our numerical simulations in which all cooling 
flows that were sufficiently heated to noticeably reduce 
the central cooling M also resulted in strongly negative ap- 
parent temperature gradients (Brighcnti & Mathews 2002; 
2003). 

However, the bubble return model clearly fails if the 
bubbles are too small. As rrib decreases, the large num- 
ber of bubbles required to transport M (particularly near 
the central heating source at r^) consumes all the avail- 
able volume, entirely displacing the cooling inflow (Condi- 
tion 21). In these circumstances the heated central region 
would grow in size and would almost certainly result in 
an extended hot thermal core and d{T)/dr < 0, which is 
not generally observed in X-rays. The circulation flows we 
consider do not specify the physical nature of the heating 
process interior to the heating radius r^. Consequently, if 
the bubbles are heated too slowly, there is some concern 
that emission from r < could result in negative appar- 
ent d{T)/dr and this would need to be considered in a 
more detailed model. We have found that the problems of 
bubble congestion near and emission within cannot 
be alleviated by increasing the radius rh at which heat- 
ing occurs. Finally, bubbles of low mass also experience 
more drag, move more slowly and, if mf, is too small, can 
fail to transport the mass arriving in the cooling inflow 
(Equation 22). 

These difficulties with small bubbles clearly raise con- 
cerns about the possible fragmentation of larger bub- 



bles. However, at the X-ray resolution of the observa- 
tions we have so far, the bubble-like disturbances have 
sizes rb/r ~ 0.05 — 0.3 that are consistent with those re- 
quired in our circulation flows. It is also noteworthy that 
the X-ray images showing bubbles and surface brightness 
fluctuations are rather similar on all scales, from elliptical 
galaxies to clusters like Perseus. The azimuthally random 
character of bubble or density fluctuations at all radii in 
the hot gas also seems to be rather universal. The az- 
imuthally random orientation of bubbles may be difficult 
to explain with numerical simulations in which energy is 
deposited by a fixed jet; the jet orientation would need 
to change rapidly due to precession or other similar black 
hole disturbances. Finally, the X-ray data so far indicates 
that the strong (or complete) reduction in the cooling rate 
M occurs universally in all hot gas flows, whereas only 
a small (but extensively observed) fraction of these flows 
contain powerful jet-like radio sources. 

Circulation flows are also limited in spatial extent, and 
this limitation becomes particularly acute as the bubble 
mass decreases. In all the galaxy-scale circulation flows 
that we considered, the circulation radius Vc is limited 
to r < 25 kpc which may be comparable to observed re- 
gions having a multi-temperature spectrum. For example 
in NGC 5044 Buote et al. (2003) find evidence of signifi- 
cant temperature fluctuations out to ^ 30 kpc. Although 
small bubbles arc disadvantageous for a variety of reasons, 
there is no obvious mechanism that favors the production 
of larger bubbles or guarantees that larger bubbles do not 
fragment. This must be viewed as an essential shortcom- 
ing of cooling flows with bubble return and the notion that 
M is small or zero because of AGN heating. While it is 
likely that the circulation radius rc can be increased be- 
yond 25 kpc if the bubbles are supported with cosmic 
rays or magnetic fields, such bubbles may be unable to 
carry the requisite M away from the heated core. 

In our circulation flows, the kinetic energy and momen- 
tum of the rising bubbles is derived ultimately from the 
gravitational field and some of this energy and momentum 
is transported to the surrounding gas. The momentum ex- 
change reduces the velocity of inflowing gas, and may be 
analogous to the entrainment of local gas in many of the 
numerical simulations of individual large bubbles reviewed 
in Section 2. Bubbles can directly heat the surrounding 
gas if they thermally dissolve into the interbubble gas be- 
fore reaching the radius Tc of equal entropy. While we 
do not consider such thermal mixing here, it would almost 
certainly heat the interbubble fiow and result in unfavor- 
able apparent temperature gradients. Some limited heat- 
ing of the interbubble gas by adiabatic bubbles - either by 
drag or PdV mechanical heating or by entropy mixing - 
may account for the enhanced non-gravitational entropy 
observed in galaxy groups (e.g. Pomnan et al. 1999). 
However, the efficiency (e^ and Cpdv) of such distributed 
heating would need to be fine-tuned to be consistent with 
the flow of both mass and entropy to large radii as required 
in circulation flow models and to satisfy observational con- 
straints on the apparent temperature profile {T)(r). 

The idealized circulation flows we describe here could 
be extended to include cooling inflow from r > Tc, con- 
stant mass loss from stars, bubble fragmentation, cosmic 
ray pressure in the bubbles, a range of heating factors h 
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for the bubbles, bubble-flow thermalizatfon, and a distri- 
bution of bubble masses nib in each flow. It seems likely 
however that the filling factor and mass flow constraints 
(Equations 21 and 22) must still be valid for any flow in 
which all or most of the mass flow M is carried outward by 
heated regions. It is likely that our less complicated circu- 
lation flows contain essential features that are generic to 
all centrally heated flows designed to greatly reduce the 
cooling rate M. 

The idealized circulation flows we describe here could 
only apply for a limited time because of secular increase 
in the gas mass as new gas enters the circulating region. 
For example, if all the gas ejected from evolving red giant 
stars merges into the hot phase, the characteristic time 
for the hot gas density to increase appreciably in NGC 
4472 is U ~ /9/Q!*p* « 4.2 (rc/lO kpc)!-^^ Gyr where 
a, = 4.7 X 10~^° sec~^ is the current specific rate of stellar 
mass loss. Moreover, in many galaxy/group scale cooling 
flows, including NGC 4472, the hot gas extends far beyond 
the circulation radii rc we calculate here. Therefore, the 
cooling inflow of gas crossing radius Vc provides another 
secular increase in the mass of gas in the circulation re- 
gion. As the mass of circulating gas slowly increases, the 



circulation flow M is also expected to increase and the 
amoimt of heating h must decrease to maintain the circu- 
lation. Ultimately, a cooling event is expected (Equation 
22). When the central M is strongly reduced in our nu- 
merical simulations of centrally heated cooling flows, the 
flow cools episodically at larger radii in the flow (Brighenti 
& Mathews 2002; 2003). It remains to be determined if 
intermittent cooling events can agree with the X-ray ob- 
servations, but the observations we have so far are consis- 
tent with no cooling at all. Finally, the hot bubbles also 
carry Type la supernovae (SNIa) iron enrichment out to 
Tc, broadening the iron abundance profile relative to the 
stars, which agrees qualitatively with observations. The 
iron from SNIa acts like tracer particles that can reveal 
the direction and magnitude of the radial flow of hot gas 
(Mathews & Brighenti 2003). 
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Fig. 1. — Panels a, b and c show respectively the gas density, temperature and flow velocity for a bubble-free cooling flow from Tc = 10 
kpc to r/i = 1 kpc. These solutions of Equations (4) and (5) are designed so that the gas density matches the observed density in NGC 4472 
at Tc and r^. The observed temperature and density profiles for NGC 4472 are shown with dot-dashed lines. 

Fig. 2. — Six circulation flows with no heating from bubble-flow interactions (e = 0). Panels a, b, c and d refer to flows with bubbles of 
mass mj, = 3 X 10^ Mq. Panels a, b, and c show respectively the gas density, temperature and velocity in the cooling inflow {solid lines) 
and the gas density, temperature and velocity of the outward moving bubbles (dashed lines). In panel a the observed density profllc in NGC 
4472 is shown with a dotted line. Panel b shows the emission- weighted local temperature (T) of the flow plus bubbles (dotted line). Panel d 
shows the volume filling factor f(r) of the cooling infiow (dashed lines) and the radius of the bubbles i't,(r) (dash-dotted lines). The pair of 
flows shown in panels a - d correspond to heating factors h = 3 and 3.6. Flows with larger h extend to larger radii r'c. The second column of 
panels (e - h) shows the same information for flows with bubbles of mass mi, = 10^ Mq with and h = 3 and 6.5. Finally, the third column of 
panels (i - I) shows the same information for flows with bubbles of mass rrn, = 10® Mq with h = 3 and 6. 

Fk;. 3. — Four circulation flows with heating from bubble-fiow interactions at various drag efflccncies (with epij„ = 0). The labeling 
scheme is the same as in Figure 2. All solutions correspond to bubbles of mass ruf, = 10^ Mq and heating factors h = 6. The four solid curves 
in panel a show cooling inflows with progressively larger heating efficiencies, e = 0, 0.1, 0.3 and 0.7. The extent of the flow Tc increases with 
the heating efficiency e^, which allows each solution curve to be identified with the appropriate efficiency. 

Fig. 4. — Three circulation flows heated by expanding bubbles with Vcixious efficcncics Cpdv 

(with ed = 0). The labeling scheme is the same 

as in Figures 2. All solutions correspond to bubbles of mass = 10^ Mq and heating factors h = 6. The three solid curves in panel o show 
cooling inflows with progressively larger heating efficiencies, Cp^v = 0.1, 0.5 and 0.7. The extent of the flow rc increases with the heating 
efficiency Cd, which allows each solution curve to be identified with the appropriate eflaciency. 
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